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The last missing correction to the muon lifetime in the Standard Model at 0(ct 2 ) coming from 
gauge and Higgs boson loops is presented. The associated contribution to the parameter Ar in the 
on-shell scheme ranges from 6 x 10 -5 to —4 x 1CF 5 for Higgs boson masses from 100 GeV to 1 TeV. 
This result translates into a shift of the W boson mass which does not exceed ±1 MeV in the same 
range and amounts in particular to approximately —0.8 MeV for a 115 GeV Higgs boson. 
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The muon decay lifetime (r M ) has long been used as 
an input parameter for high precision predictions of the 
Standard Model (SM). It allows for an indirect determi- 
nation of the mass of the W boson (Myy), which suffers 
currently from a large experimental error of 39 MeV , 
one order of magnitude worse than that of the Z bo- 
son mass (Mz). A reduction of this error by the Large 
Hadron Collider (LHC) to 15 MeV @ and by a future 
linear collider to 6 MeV || would provide a stringent test 
of the SM by confronting the theoretical prediction with 
the experimental value. 

The extraction of Mw with an accuracy matching that 
of next experiments, i.e. at the level of a few MeV ne- 
cessitates radiative corrections beyond one loop order. 
Large two loop contributions from fermionic loops have 
been calculated in jl . The current prediction is affected 
by two types of uncertainties. First, apart from the still 
unknown Higgs boson mass, two input parameters in- 
troduce large errors. The current knowledge of the top 
quark mass results in an error of about 30 MeV S , which 
should be reduced by LHC to 10 MeV and by a lin- 
ear collider even down to 1.2 MeV. The inaccuracy of 
the knowledge of the running of the fine structure con- 
stant up to the Mz scale, Aa(Mz), introduces a further 
6.5 MeV error. Second, several higher order corrections 
are unknown, of which only one at the 0(a 2 ) order. This 
lacking contribution comes from diagrams with no closed 
fcrmion loops. It has been previously estimated to be of 
the order of the square of the one loop bosonic correction 
0. This, however, is unjustified since both have a dif- 
ferent dependence on the Higgs boson mass (logarithmic 
vs. quadratic). In this Letter, the question of the exact 
size of this contribution is finally settled. 

The muon decay is naturally described in the language 
of effective field theory. The process' dynamics are driven 
by the leptons, the five light quarks, and the two mass- 
less gauge bosons. Even then the momentum scale set by 
the muon mass leads to a strong decoupling of the r lep- 
ton and of the heavier of the five quarks [[| . The heavy 
particles W, Z, the Higgs boson and the top quark gen- 



erate point interactions. A general effective lagrangian 
assumes the form 

£ C ff = £qed + £ Q 5 CD + Yl (m^)" "' ^ 

(5) 

where Cqed and Cq CD are the bare QED and the five 
flavour QCD lagrangians, % n are composite operators of 
the light fields and C % n are the respective dimensionless 
matching coefficients, obtained by comparing the Green 
functions of the full theory with those of the effective one. 
The mass of the W boson has been chosen as the heavy 
scale. 

Due to the left-handed nature of the charged current 
in the SM, the only operator of dimension six relevant to 
muon decay is the four-fermion interaction of the Fermi 
model 

O f = e 7 Q (l- 7s)M® VfoCl -7sH, ( 2 ) 

with /i, e, v e denoting the spinors of the muon, elec- 
tron and their neutrinos. Here the so called charge- 
conserving form has been used (Fierz transformed SM 
amplitude). The respective matching coefficient is tradi- 
tionally parametrised as 

C F = M^ = ^(1 + Ar), (3) 

where Gf is known as the Fermi constant, whereas Ar 
represents the higher order corrections and vanishes at 
the Born level. The one loop corrections in the effective 
theory given by Eqs. [j], || and || have been calculated in 
and S, while the two loop in |J. 

The framework for construction of effective lagrangians 
together with a proof has been given in pi . It is based 
on expansions of individual Feynman diagrams. The 
method is presented on the example of a specific graph 
contributing to bosonic corrections to muon decay at two 
loop order in Fig. |l|. The terms in the large mass expan- 
sion can be classified according to the scale of the loop 
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FIG. 1: A typical muon decay diagram a) and the contribu- 
tions to its large mass expansion according to the momenta 
b) fci-soft, fe-soft; c) soft-hard; d) hard-soft; e) hard-hard. 

momenta. When both momenta are "soft" (<C Mw), as 
in b), the propagators of the W and Z bosons are ex- 
panded leading to a correction of order a/M^ in the ef- 
fective theory. For one momentum "soft" and one "hard" 
(~ Mw), as in c) and d), corrections of either order, 
a/M^y or l/M^r in the effective theory, are generated. 
The contribution to the matching coefficient comes only 
from the region where both momenta are hard, as in e). 
In this case all of the light particle masses and momenta 
should be put to zero. By these arguments it can be 
shown that Ar can be obtained by simply taking the sum 
of all the diagrams and putting all external momenta and 
light masses to zero. The procedure should generate no 
spurious infrared divergences, while the physical diver- 
gences connected with the photon should be contained 
in the corrections of the effective theory. As known the 
Fermi theory corrections are finite, therefore the Ar cor- 
rection obtained as above should also be finite. 

Previous calculations of Ar have been based on a dif- 
ferent method of factorisation originally devised in [jll| . 
This procedure consists of subtracting from the infrared 
divergent SM diagrams the respective Fermi theory di- 
agrams in Pauli-Villars regularisation. The difference is 
well defined in the limit of zero light masses and external 
momenta. It turns out, however, that the QED Ward 
identity, which is responsible for the finiteness of the cor- 
rections in the Fermi theory, implies in this case the van- 
ishing of the sum of the subtracted diagrams. This proves 
that both procedures are equivalent. 

The evaluation of two loop corrections to a four- 
fcrmion process requires the full second order renormal- 
isation of the SM lagrangian in all but the Higgs sector, 
where first order suffices. The comparison with experi- 
ment imposes the use of on-shell parameters for the final 
result. Throughout this work the on-shell scheme was 



used, with a procedure similar to the one described in 
The only substantial difference concerns the treat- 
ment of tadpoles. 

It is known that gauge invariance of mass counter- 
terms requires inclusion of tadpoles [[l2[ [ll| (a t the two 
loop level this has been explicitely shown in In this 

case, however, one cannot use one-particle-irreducible 
(1PI) Green functions. In order to have gauge invariant 
counter-terms and 1PI Green functions only, a special 
procedure was designed. An additional renormalisation 
constant for the bare vacuum expectation value vq, de- 
noted Z v , has been introduced and explicitely split from 
the bare masses 

v — > v Zl'\ (4) 
M^ z — M^ Z Z\I\ (5) 

The term linear in the Higgs field H in the lagrangian 

T H° = ^A(M^) 2 Zl/ 2 (Z v -l)H , (6) 
eo 

is then used to determine Z v , through the requirement 
that tadpoles are canceled. It can be proved |l2[ [ll| that 
the bare masses are gauge invariant in this case (an equiv- 
alent procedure which makes use of the effective potential 
has been used in Q). 

The calculation of the two loop bosonic contributions 
to muon decay was performed by means of a completely 
automated system. The diagram generation stage was 
done by the C++ library DiaGen |l7). The tensor re- 
duction of two loop propagator diagrams was accom- 
plished with the algorithm described in [lq| , whereas 
vacuum diagrams were treated with integration by parts 
identities |l9fl . For algebraic manipulations, the program 
FORM [ pCf was used. The two loop two-point inte- 
grals were numerically evaluated with single integral rep- 
resentations of the package S2LSE The latter was 
modified for quadruple precision, which was needed due 
to large cancellations (independent terms grow as Mjj, 
while the result behaves as Mfj). 

The size of the software required several tests. The 
following algebraic checks were performed 

• ultraviolet and infrared finiteness, by cancellation 
of poles in dimensional regularisation, 

• gauge invariance, by independence of the three 
gauge parameters of the general gauge for the 
SM, 

• Slavnov- Taylor identities for two-point functions, 
as given in jl8), both for on-shell integrals and by 
expansion in the external momentum to second or- 
der. 

Several numerical tests were also done 
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FIG. 2: Ar as function of the Higgs boson mass (solid line) 
and the leading term in the large Higgs mass expansion 
(dashed line). 



FIG. 3: W boson mass shift generated by two loop bosonic 
corrections. 



• all of the master integrals were evaluated indepen- 
dently by means of deep mass difference and large 
mass expansions, 

• each of the two-point on-shell diagrams was calcu- 
lated separately with the help of small-momentum 
and different large-mass expansions, 

• the result of (IJ] for the W and Z mass counter- 
terms was reproduced to precision dictated by the 
order of the expansions contained therein. 

A detailed description of the methods used is relegated 
to a subsequent publication Jl5| ]. 

The variation of the bosonic contribution with Mh 
subtracted at Mjj = 100 GeV has also been compared 
with the one presented in Q . A discrepancy which grows 
up to ~ 16% for Mh = 1 TeV, has been found. 

All the numerical values have been obtained for the 
following coupling and mass parameters Q 



(7) 



a' 1 = 137.03599976(50), 
G F = 1.16639(1) x 10~ 5 GeV" 2 , 
M w = 80.423(39) GeV, 
M z = 91.1876(21) GeV. 

(21 

The two loop bosonic correction to Ar, denoted Ar y bo ' s , 
is shown in Fig. ^ as function of the Higgs boson mass 
Mh- Its behaviour is very smooth apart from the range 
between 150 GeV and 190 GeV, where two peaks are 
located. Their positions at ~ 161 GeV and ~ 182 GeV 
correspond to the W and Z pair production thresholds 
introduced through the on-shell definition of the Higgs 
boson mass. 

It is interesting to study the correction for large values 
of Mh- In fact, two different leading M H terms have been 
given in p2| and p3|, both based on the same results 



of p4| . In this last work a renormalisation scheme has 
been used in which M H terms occur only in two-point 
functions. The leading behaviour can then be obtained 
from the formula 



(2) 
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(8) 

where Se is the charge renormalisation and AMw,z rep- 
resent the shift of the W and Z self energies from zero 
external momentum up to the on-shell value. Ap has 
been introduced in |2f|. With the help of p4j , the result 
reads 



Ar 
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8M W 



(9) 
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where CI2 is the Clausen function and Cl2(7r/3) = 
1.0149416. Eq. ^ is in agreement with j22), whereas in 
p3[ the charge renormalisation has not been included. 
The leading Higgs mass dependence is depicted in Fig. |^. 
The difference between the full result and Eq. ^ for values 
below 1 TeV shows the importance of subleading correc- 
tions. It has been checked that for Mh > 5 TeV both 
functions agree very well. 

With the help of Ar the W boson mass can be esti- 
mated by inverting Eq. || 



M w = M z 



\ 




V2G F M; 



-(1 + Ar). (10) 



This equation should be solved recursively, since Ar also 
depends on Mw- In the case of small corrections one 
can, however, obtain the additional mass shift by means 
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of an expansion. For parameters' values given in Eq. [7], 
a second order Taylor series yields 

AM W = -(1.491 + 1.779 Ar) x 10 4 Ar^ [MeV], (11) 

where Ar represents the dominating correction, which 
is composed of one loop electroweak, two loop fermionic 
and O(aas) QCD contributions. Since Ar does not ex- 
ceed 5% ||, the third order in the Taylor series would 
contribute less than a percent which is comparable to 

(2) 

the inaccuracy induced by the dependence of Ar£ ^ on 
the W mass. 

The W mass shift, Eq. with Ar neglected, is plotted 
in Fig. |3|. The dropped term implies an error of less 
than 6%, which is, however, completely negligible in view 
of the size of the total effect, which does never exceed 
±1 MeV. 

In conclusion, the complete 0(a 2 ) bosonic contribu- 
tions to muon decay have been calculated in the on-shell 
scheme, and shown to lead to a small W boson mass shift 
below 1 MeV. A framework for evaluation of Ar which 
does not explicitely refer to the QED corrections in the 
Fermi theory has been presented and its equivalence to 
the method based on subtraction of Pauli-Villars regu- 
lated diagrams demonstrated. A procedure of obtaining 
gauge invariant mass counter-terms without explicitely 
including tadpole graphs has been developed. 
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